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We review the exact results on the various critical regimes of the antiferromagnetic te- 
state Potts model. We focus on the Bothe Ansatz approach for generic Q, and describe 
in each case the effective degrees of freedom appearing in the continuum limit. 
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1. Introduction 



The Q-state Potts model-'^, defined in 1952, is a lattice model for classical mag- 
netism. Since the 1970's, it has been used as a laboratory to develop theoretical 
methods of Statistical Mechanics, such as Kramers- Wan nier duality, Yang-Baxter 
integrability, Bethe AnsatJ^I^ Coulomb-gas formalisnJ^El^ Conformal Field TheorjEl 
(CFT). Most studies on the Potts model deal with the ferromagnetic critical point, 
described by a Coulomb-Gas CFT or minimal CFT at the rational values. In con- 
trast, the present review is concerned with the antiferromagnetic (AF) regime of 
the Potts model, which also contains exactly solvable critical points. Although it 
is based on the same simple lattice model, this regime exhibits interesting physical 
features. Indeed, non-unitarity of the equivalent vertex model allows a CFT with 
central charge c > 1, and non-compact degrees of freedom. 

Let us introduce the subject with a short historical notice. In 1982, Baxter first 
discovered^ the location of the AF/paramagnetic transition by a mapping to a 
staggered integrable six-vertex (6V) model, and calculated the free energy density. 
Later on, Saleur describecffi the complete phase diagram for < Q < 4, emphasing 
the role of the Beraha numbers Q — 4 cos^ j for integer t. At these values, he 
gav the CFT description of the AF/paramagnetic transition, relating it to the 
Zfc parafermion This transition corresponds to a staggered integrable model, 

but it appearetJI^ as a solution of the basic Yang-Baxter Equations (YBE) for the 
{Na,Np) Potts model. A comple te CF T description of this transition for both 
generic and Beraha Q was achieved ^^ l ^^ l, especially through the study of the Bethe 
Ansatz equations; it consists in one compact boson (f>i with Q-dependent radius 
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and one non-compact boson (f>2- Besides the isotropic critical regimes mentione d so 
far, a mixed ferro-antiferromagnetic anisotropic regime was defined and solved^. 
Its continuum limit consists in two compact bosons: 0i with a Q-dependent radius, 
and 4>2 with a fixed radius, such that (j)2 decomposes into two Ising models. 

This review is organised as follows. Section [2] is a reminder on the Potts model 
and its phase diagram for generic < Q < 4, and describes some related problems 
relevant to Statistical Mechanics, but also Condensed Matter Theory. There are 
three critical regimes for the Potts model with AF interactions. In Section [3l we 
summarize the methods used to study the critical points of the Potts model for 
generic values of Q. Sections IH [5l |6] then apply these methods respectively to the 
BK phase, the AF transition line and the anisotropic critical regime. Section [7] 
concludes by mentioning some aspects that are not treated in detail by this review, 
together with some interesting open problems. 

2. The Potts model and related Physics problems 
2.1. Definitions and integrability properties 

Let £ be a planar lattice. The Q-state Potts model on C consists of spins Sj living 
on the vertices of £, and taking the values {1,2,..., Q}. The Boltzmann weight for 
a spin configuration is 

W[{S,}]=expljJ2Ss.,sA , (1) 
V (^i> / 

where 6 stands for the Kronecker symbol and the product is on all edges of £. A 
positive (resp. negative) value of J favors configurations where neighboring sites 
have equal (resp. distinct) spins. This defines the ferromagnetic and AF regimes. 

The graphical expansion of the partition sum gives the Fortuin-Kasteleyn (FK) 
cluster modeP^ (see Fig[T]): 

^Pott. = Q''"^^"'*''^ ■■= Zfk{C, Q,v), v.^e' -I, (2) 

where C{G) is the number of connected components (clusters) of the subgraph G, 
and £(G) is the number of edges in G. In this representation, Q plays the role of 
a real parameter, and we are mainly concerned with the regime < Q < 4. Using 
Euler's relation, the partition sum can be expressed as 

Zfk(£,Q,w) = const X ^ Q^^^M-I = const X ZFK(>C',Q,g/w), (3) 

G'C£' ^ " ^ 

where C is the lattice dual to C. We refer to the transformation ^ as Kramers- 
Wannier duality. 

The FK model can be, in turn, mapped to a loop model based on the Temperley- 
Lieb (TL) algebrePI (see Fig [1]) . The TL loop model lives on the medial lattice 



I 



I 
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Fig. 1. Example FK cluster configuration on the square lattice C, and corresponding loop con- 
figuration. Direct (rcsp. dual) clusters arc shown as thick full (resp. dotted) lines, and clusters 
consisting of a single site as dots (resp. crosses). 



consisting of the midpoints of the original lattice jC. A cluster configin-ation defines 
imiqucly a loop configuration, by asking loops to separate direct and dual clusters. 
We have the relation 



(4) 



I n = 

loop config. \x = vl^ 

where c{H) is the number of closed loops in H, and X{H) is the number of medial 
sites where the loops do not cross the original edge of C. Although the Potts model 
can be defined on any planar lattice, we restrict the subsequent discussion to the 
square lattice for simplicity. 

Consider the Potts model on a strip on width L sites (with L even), with coupling 
constant Ji for horizontal edges, J2 for vertical edges, and corresponding weights 
X1.X2 in the loop model. The transfer matrix for loops can be built with the TL 
algebra TLi(n). This algebra is generated by the operators Cj for j — 1, . . . , L — 1, 
satisfying: 



for \j — fc| > 1. 

The transfer matrix parallel to one axis of the original lattice £ then reads 



n(i+ 



Xl Bj 



J even 



j odd 



(5) 



(6) 



Our regime of interest is < n < 2. We introduce the crossing parameter 7 and 
the parameter t, defined by 



V^=2cos7, 0<7<J, t:=-. 

2 7 



(7) 



The algebra TLi(n) provides a solution to the YBE, with one spectral parameter 



u: 



Rjj+i{u) = sin(7 — u) 1 + sin(u) Cj . 



(8) 



I 



I 
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Table 1. The critical regimes of the Potts model on the square lattice. The last column 
gives the generic value of the central charge for the loop model. 



self-dual 


< n < 7 


XI, X2 > 


ferromagnetic critical 


c- 1 « 




7 < M < 7r/2 


xi,X2 < 


BK phase 




staggered 


7 < M < tt/2 
< M < 7 


xi,X2 < 
xi > 0, a;2 < 


AF critical 
anisotropic critical 


^ _ o 12 

t{t-2) 



2.2. Phase diagram 

Let us recall the locatiorPEl of the critical points of the square-lattice Potts model. 

• From ^ and ((H]), we see that the self-dual line xiX2 = 1 corresponds to a 
homogeneous loop model with 

sinu 1 

Xl = r , 2^2 = ■ (9) 

sin(7 — u) xi 



• The only other integrable case respecting the parit y of site s in ([6]) is when 
spectral parameters alternate between u and u + H6 | 12 | 13 1 



sinw cos(7 — u) 

XI = —, r , X2 = . 10 

sm(7 — u) cosu 

Each of these two cases contains two critical regimes, corresponding to different 
values of the spectral parameters: see Table [1] The three isotropic critical regimes 
are shown in Fig. [21 



2.3. Related problems of Statistical Mechanics and Condensed 
Matter Theory 

Using the FK formulation, one can relate the Potts model to various statistical 
models through exact equivalences. 

• Percolation. At Q = 1 , the FK model reduces to bond percolation on the 
lattice £. 

• Uniform spanning trees. In the limit {Q ^ 0,v ^ VQ), using Euler's 
relation, 

Q-i(^+^)Zp^{£^ Q, y/Q) — > # spanning trees on £ , (11) 

where J\f is the total number of sites in £. So the continuum limit of the 
uniform spanning tree problem^ is determined by the ferromagnetic Potts 
transition at Q — > 0. 



I 



I 
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Fig. 2. Phase diagram of the isotropic Potts model on the square lattice. The self-dual line v = 
is shown as a thin line, and the AF critical line d = — 2 it y/A — Q as a thick line. The BK 
critical phase, governed by the self-dual branch v = —\/Q, is enclosed in the thick line. Arrows 
represent the renormalization-group flows in the vicinity of fixed points. 



• Spanning forests. In the vicinity of Q = 0, one can define another hmit: 
(Q — > 0,w = Q/w), where w is constant. Using the same arguments as 
above, we get the hmit: 



iQ/w)-<^^+''>ZM^,Q,Q/w) V (12) 

Q-i-O ^ — ' 



where T{C) is the set of spanning forestall on £, i.e. spanning subgraphs 
with no internal face. The parameter w thus plays the role of a tree fugacity 
in the spa nning forest problem. The phase diagram of the Potts model 
predicts^ the behavior of spanning forest as a function of ?i;: if uj > or 
w < —4, the model is not critical (paramagnetic phase); if —4 < u; < 0, 
the model is described by a Coulomb-Gas theory (BK phase); exac tly at 
w = —4, the model is described by an OSP(2|2) supersphere a-model^. 
• Restricted Solid-On-Solid (RSOS) models. At the Beraha numbers 
-JQ = 2 cos J with integer t > 4, the TL generators can be represented as 
local operators in the At-i RSOS height The latter consists 

of local variables {hj,j — 1, . . . , L} subject to the local constraints 



/i, e {i,...,i-i}, \h,-hj+i\ = i. (13) 
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The operators e,- are defined as 



sm sm ^ 



ej\hi...hL) =Sh^_,,h,+, J2 . ih,+, * \hi...h' ...Hl) . (14) 

The RSOS model may describe the statistics of a two-dimensional interface 
in three-dimensional space. 

A second approach is to take the very anisotropic limit m — >■ 0, where the 
transfer matrislfj generates a local one-dimensional quantum Hamiltonian based on 
the TLi(n) algebr£0, and consider the spin-i representation. The self-dual © and 
staggered (fTO|) transfer matrices respectively give the HamiltonianPIEl 

L 

iJ.d=±^e,, (15) 

L 

iJstag = ± ^(-n Cj -t- ejcj+i + ej+iej) . (16) 
i=i 

The global minus sign in H^d (resp. i?stag) corresponds to the ferromagnetic (resp. 
AF) critical point; the plus sign to the BK phase (resp. the anisotropic critical 
point). 

The TL generators can be represented as rescaled projectors in the Ug(Sl2) 
quantum algebra"'^^, where q = e*"*". Consider a collection of L spin-i variables 
{cTj}. In terms of the Pauli matrices ctJ'^'^, the self-dual Hamiltonian describes an 
XXZ quantum spin chairPl 



L 

I . 



1 ^ 

If we denote Pf^j+i the projector of {oj -\- (Jj+i) onto total spin zero, then the 
operators Cj := "^-Pj "+1 satisfy the TLi(n) relations ([5]). Furthermore, the projector 
of [aj + (Tj+i + CTj+2) onto total spin | is given by 

p(3/2) _ , (cj-ej+i + Cj+iej) - 7i(ej + e^+i) 
- + n2 - 1 

In terms of the Ug(Sl2) projectors, the staggered Hamiltonian readi^ 

7j.tag = ± E + - • (19) 

If the plus sign is chosen and n > 1 , the second term in ([TOl) favors a totally dimer- 
ized state (a product of L/2 Uq(Sl2) singlets), while the first term gives additional 

^ We consider here the transfer matrix parallel to one axis of the medial lattice M. 
^ We use a periodic variant of T'Li^{n), with ej±ij := Sj. 



I 
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energy to dimerized pairs of neighbours: the spin chain is thus subject to frustra- 
tion. So -ffstag is an integrable, g-deformed generahzation of the weh-studied SU(2) 
frustrated spin 



3. General approach to the solution 

3.1. Equivalent vertex model, twisted BC 

We consider the Potts model on the lattice L, embedded in a cylinder. Its exact 
solution is obtained through the sequence of mappingj^ 

(Potts on C) — > (FK on C) — > (TL on M) — > (6V on M). 

The resulting 6V model is homogeneous for self-dual Potts, and staggered for the AF 
and anisotropic critical lines. The last mapping is done by orienting independently 
each loop, and giving it a weight e^*''' according to its orientation, so that the total 
weight for a loop is n = 2 cos 7. The complex weight e^'^'^ is distributed locally 
in the 6V model, by counting the local turns of each loop fragmentP. However, 
through this mapping, the loops which wind around the cylinder {non-contractible 
loops) recieve an incorrect weight n = 2, because their total rotation angle is zero 
in both orientations. To correct this, we consider the 6V with a seam along the axis 
of the cylinder: an arrow which crosses the seam from left to right (resp. from right 
to left) gets a factor e*"^ (resp. e"*'''). Non-contractible loops then get a weight 

n — 2 cos if . (21) 

The Potts ground state corresponds to = 7, but it can be useful to consider other 
values of (f, e.g., to compute the magnetic exponent of the Potts model (see below). 



3.2. Analysis of the BAE 

The integrable 6V is solvable by Bethe Ansatz. In all the cases we consider here, 
the Bethe roots sit on horizontal line(s) in the complex plane, and are described by 
a root density p in the large-L limit. One may extract from the BAE two types of 
data which are relevant to the associated critical theory: 

• The physical properties of elementary excitations (particle-hole excitations 
in the root distribution) are derived easily by Fourier transform. This leads 
to the scattering theory of these quasi-particles both at criticality and in 
the vicinity of the critical point. 

• The finite-size corrections to the energies allow us to relate the lattice model 
to GET, where the energies and momenta are expected to behave as 

live 

EQ^e.^L~ — , (22) 
bL 

'2iTTU — 

E,^-f^c^Eo + —{h + h), (23) 
khrh = Y^h-h) mod 2^, (24) 
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where Eq is the ground-state energy, is the energy density, v is the 
Fermi velocity, c is the central charge and E/^ hT^hh the energy and 
total momentum associated to the primary state with dimensions {h,h). 
The computation of the finite-size corrections involves the Wiener-Hopf 
technique, described in this context by Yang and Yang^. 



3.3. Numerical transfer-matrix computations 

Analytical results may be checked numerically by computing the dominant eigenval- 
ues of the transfer matrix in various sectors or, equivalently, the lowest eigenvalues 
of the Hamiltonian. Through ([22|) - ([23l) . the eigenvalues provide numerical estimates 
for the central charge and conformal dimensions. 

Using the TL formulation, one usually computes systems of up to i = 20 sites, 
and the precision obtained for the exponents is about 10^^ (except in the presence 
of logarithmic corrections, typically at 7 — ^ 0). 



4. The Berker-Kadanoff phase 

4.1. Relation to the XXZ spin chain 

The two branches of the self-dual line (|9|) are related to each other by {xi,X2) — >■ 
(— .Ti,— 2:2), or equivalently by y/Q — > — \/Q (see (O). We denote /i := tt — 7, so 
that —\fQ = 2cos/i. In the very anisotropic limit 11 — > 0, the BK critical line is 
thus equivalent to an XXZ spin chain with twisted BC: 

1 ^ 

i=i 

where 

A = — cos /i , ip = ^ . (26) 

In the present Section, we will then recall Bethe-Ansatz results on XXZ, and their 
application to the study of the BK phase for generic Q. 

On each site of the XXZ chain, we consider an up spin as an empty site, and a 
down spin as a particle. The Bethe Ansatz Equations and eigenvalues for r particles 
read: 

" sinh(f ~ A,) 
sinh(f + Aj) 



(25) 



L 



sinh(z/j, — Xj -\- \i) 
sinh(i/x -I- Xj — Xi) 



-2 sin^ /i 



cosh Xj — cos /X 



(27) 
(28) 



I 



I 
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4.2. The ground state and the quasi-particle picture 

We take the logarithmic form of the BAE: 

r 

Lk{Xj) = 2ttIj + 2<^ - ^ e(Aj - A^) , (29) 

£=1 

where 

u\\ -1 sinh(ii^ - A) sm\\{iii + A) 

fc(A) :=-»log ^.^^^^^^^^^ , 9(A) :=-^log ^.^^^^^^_^^ , (30) 

and the "Bethe integers" Ij are half-odd integers (resp. integers) if r is even (resp. 
odd). In the large-L limit, the Bethe roots \j are described by a root density /o(A), 
satisfying the linear integral equation 

/•A+ 

fc'(A) = 27rp(A) - / dv p{v)K{\-u), (31) 



where the kernel K is given by K 0', and A± are the extremal values of the 
A/s. 

The ground state for (/? = corresponds to a symmetric root distribution with 
r = L/2 and A± = ±cx), so (pij) becomes solvable by Fourier transform. We use the 
notations 

f + 00 



f{Lo):= dA/(A)e-\ if * 9){X) := f{v)g{X ^ v) . (32) 

J —oo J —oo 

The ground-state solution is 

Any physical quantity is given by an integral of the form 

Aip) = dA p(A)a(A) , (34) 

where, e.g., a = 1 for the particle density, a = k for the total momentum, etc. A 
hole at position A^ in the root distribution produces a variation in A, given by the 
dressed quantitji^a: 



A{p) - A{p,) = 



1 -ka 

2tt 



{\h):=-^a(Xh)- (35) 



In particular, the dispersion relation of the quasi-particles is given by: 

r , , ^ . / , ttA \ , , TT sin u 

k{X) = 2Arctan tanh — , e(A) = ^-3- . (36) 

V 2/i/ ^cosh^ 

Close to the Fermi levels kp = ±-|, this becomes a linear relation: 

e ~ v\k — kp] , v = —. (37) 



I 



I 
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To complete the S'-matrix picture, one defines the dressed scattering amplitudes as 
follows. Consider a root distribution with particles of density p and holes of density 
Ph- The BAE then read 



k' = 2tt{p + Ph) - K * p . 
Introducing $ (27r - if ★ {-2t:K), we express the BAE as 

k' = 2Tr{p + ph) - 



(38) 



(39) 



In the above form, the BAE describe the scattering of holes with momentum k and 
scattering amplitude 



S{\, A') = exp 



A-A' 



exp 



dv <J>(z/) 

° du sinh (^p — ^'^ uj sin(A — X')u 



2 cosh ' 



sinh ^^-^^uj 



(40) 

(41) 



4.3. CFT spectrum 

In this paragraph, we use the BAE to derive the CFT associated to the untwisted 
{if = 0) spin chain. In the next paragraph, we will then reintroduce the twist to 
compute correctly the critical exponents in the Potts model. 

As explained above, the CFT spectrum is obtained through the finite-size cor- 
rections to the energies. The central charge is give n b y the ground-state energy 
Eq. In this case, the corrections are easily computed^ll using the Euler-Maclaurin 
formula. This gives 

r+°° TTV 

Eo^L dX po(A)e(A) - — , (42) 
J-oo 6L 

and so the central charge is c = 1. 

The states corresponding to primary operators in the CFT are combinations 
of "magnetic" and "electric" excitations. A magnetic excitation of charge m G 
Z consists in removing m roots from the ground state, and arranging the Bethe 
integers Ij symmetrically around zero. An electric excitation of charge e e Z is a 
global shift of all Bethe integers: Ij Ij — e. The associated energy corrections 
may be computed by the Wiener-Hopf technique'^. This results in the expression 

On the other hand, the total momentum for the (e, m) excitation is 

2tt ^ 27rem , „ 

ke,m = ^ = — ^ 27r . (44) 



I 
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From i?e,m and fce,m, we deduce the conformal dimensions: 

1/e \^ - l/e \^ TT — ^ 

= — 1- m.. n h„ .„ = — m.. /o _ o = 



/le,m = ^ + , /le.m = ^ [^-^ - , g = — — . (45) 

We call the corresponding primary operator in the CFT. One can also identify 
the Bethe root configurations associated to the descendants of Oe,m- the action of 
the Virasoro generator L_„ (resp. Z_„) amounts to shifting the highest (resp. 
lowest) Bethe integer by n (resp. —n). 

The central charge c — 1, the primary dimensions ()45p and the Virosoro descen- 
dants obtained by this analysis of the BAE correspond to a compactificd boson (or 
Coulomb gas) CFT with action^^ 

j A^x {\/(t)f , (t) = (j) + 2TT. (46) 

In this theory, a magnetic excitation is a dislocation defect of amplitude 27rm, 
whereas an electric excitation corresponds to the insertion of the vertex operator 
exp(ze</)). 

The continuum partition sum on a torus of modular ratio t has the form: 
where t] is the Dedekind function. 



4.4. Effective central charge and critical exponents 



In the presence of a twist Lp, the conformal dimensions (j45j) are modified by the 
change of electric charges 

e e + - . (48) 

TT 

The effective central charge and critical exponents for the Potts model in the BK 
phase are obtained by considering the appropriate sectors for m and in the 6V 
model. 

The Potts ground state has charges (cq — ii/ir^m — 0), and hence the effective 
central charge for the Potts model is 

9 t 

The magnetic, £-\eg watermelon and thermal exponents Xh, Xi and Xt are 
defined through the relations 



Gff (r, r') := (^Ss^^s,, ^ ^) - ^ ^'^ 



l\-2X„ 



(50) 



GpXr, r') := P [r, r' are connected by i strands] ~ |r - r'| ' , (51) 
i^\v-vA-\ Xt:^2^-, (52) 

V 
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where ^ is the correlation length and Vc — —\/Q is the RG fixed point for the BK 
phase. 

• The magnetic correlation function may be written as 

G^(r,/) = (l-Q-i)^^^, (53) 

where ZH(r,r') is the partition sum of FK configurations where r and r' 
belong to the same cluster. In the cylinder geometry, we put r and r' at the 
ends of the cylinder, and the constraint for Zh is equivalent to forbidding 
any non-contractible loop: this is easily implemented by setting (p — n/2, 
so that n = 0. Hence, the magnetic exponent is 

" 2g 2g 2g ^ ' 

• The watermelon exponents correspond to higher magnetization sectors. 
Note that i has to be even when L is even. Exponent corresponds 
to the charges (e = 0, m = £/2), and hence: 

8 25 ^ ^ 

• The thermal exponent Xt corresponds to an irrelevant perturbation around 
V = —\fQ (see Fig. [2]), and thus we expect Xt > 2. In the sector m = 0, 
non-contractible loops are allowed, and the twist must be of the form if = 
7r(eo +p), with integer p. This leads to the series of exponents 

Xrip)^^-^^±f^. (56) 
The two lowest exponents in this series correspond to the values p — —2,1: 

Xt^2, = a _ 1 > 2 (57) 

5. The antiferromagnetic critical line 

5.1. The Z2 staggered vertex model 

The antiferromagnetic critical line is mapped to the staggered 6V model with an 
alternation of spectral parameters, as shown in Fig. |3l We consider this model on 
a cylinder of circumference L = 2N sites. Since it is built on the TL algebra, the 
transfer matrix has Ug(Sl2) invariance. Moreover, the YBE imply an additional 
symmetry, expressed by the "Z2 conjugation operator" C: 

C:^f[(-^]. (58) 
fj,\ cos J J 

The operator C satisfies = 1 and commutes with the transfer matrix. As a 
consequence, the eigenstates may be labelled according to this Z2 symmetry. 



I 



I 



13, 2010 1:10 WSPC/INSTRUCTION FILE paf-rev-v2 



The antiferromagnetic Potts model 13 











































1 


\ / 


\ 1 


\ 1 


\ 



f I 



Fig. 3. Spectral parameters of the Z2 staggered vertex model. 



5.2. Bethe Ansatz Equations 

The BAE and energies for the Z2 staggered model read, for even N: 

sinh(i7 + aj)]^ ^ ^_^yr-l^2^v TT sinh^(2z7 + aj--a^) ^^^^ 

fj-^ sinh i(2i7 - aj + ag) ' 



sinh(i7 — aj) 



E^± 'r^'^ ^ . (60) 

^-^ cosh 2a,- — cos 27 
j=i J 

The Bethe roots are individually defined modulo 2«7r, but the BAE are invariant 
under a global shift aj — >■ oj + iir. This refiects the Z2 invariance discussed above, 
and in fact, one can prove^that the conjugation operator C acts as follows on the 
Bethe eigenstates {ipiai, . . . , ar)): 

C|V'(ai, • ■ • , ttr)) = const x |'0(ai + jtt, . . . , + irr)) . (61) 

The expression ((60)) for the energy suggests that Bethe roots aj lie on the lines 
Ini a = ±2j, so that their contribution to the energy is negative. Hence, we divide 
the roots into two sets 

«° = A° + f, a]-A]-f. (62) 

We denote r° the number of roots A°. In terms of the A°, the BAE and energies 
read 

exp [tNk{X^)] = (-l)'-°-ie2'^ 11 ^'^P [-^©''"'(A ■ - A,^)] , (63) 

6=0,1 e=i 

a=0,lj = l a=0,lj = l 3 ' 



where: 



-1 cosh(z7 + A) 

k{X) := -I log — , (65 

cosh(i7 — A) 
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f^^tw -1 sinh(z7 - I) ^±i.,x cosh(z7- I) 

Q (A) := -»log . ■ ,. , , 9=^ (A) := -ilog ^ ■ (66) 

sinh(z7 + ^) cosn(i7 + ^) 

In these notations, the Z2 conjugation (j6T|) exchanges the indices a = 0, 1. 



5.3. XXZ states 

If we denote /i = tt — 27, then we have the identities 

sinh(f - A) 



k{\) = —i log ■ 



sinh(f + A) 



(eo + e±^)(A)^-.iog^!^^ifc^ 



sinh(i/i — A) 
-2 sin^ 



(67) 
(68) 



e(A) = — — ^ . (69) 

cosh 2A — cos 

Hence, if the sets {A°} and {Aj} are identical, the BAE and energies are equivalent 
to those of the XXZ spin chain with parameter ^ and twist Lp. 

5.4. Ground state and quasi-particles 

The discussion of Section |4] on the large-L limit of the BAE can be adapted to the 
Z2 staggered model, except we now have two coupled sets of BAE satisfying Z2 
symmetry. The logarithmic form of the BAE is 

iVfc(A°) = 2nl^ + 2^ - ^ e°(A;; " A°) - ^ e-^A," - A^) , (70) 
£=1 i=i 

Nk{X]) = 2nlj + 2^ - ^ e\X] - A°) - ^ e°(A] - A^) , (71) 

e=i i=i 

where the are half-odd integers (resp. integers) if is even (resp. odd). In the 
large- hmit, the roots lie on the intervals [A''_,A^] and [A^,A^], and the BAE 
become 

fc'(A) = 2^p'^(A) - V \ ^ Av p^{u)K^-\\ - v) , (72) 

where :— (0")'. The ground state is a particular case of the XXZ states de- 
scribed in Section [5.31 ^-nd hence the ground-state root densities are 

/(A)...(A)...,A).^^, (73) 

A physical quantity reads 

Aip°,p^)^ [ ^dA p*'(A)a(A) + / ""dA pi(A)a(A) (74) 
Ja° Ja^ 



I 



I 
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We introduce the linear combinations K'^ ± . The presence of a hole 

in the root distribution {A°} results in the change: 

-1 



^(pO,pi)-A(po,Po) = -^ 



K 

1 I * a 

27r 



(A?J:=S(A°) 



(75) 



So, from ()67p - (|69p . the elementary excitations have the same dispersion rela- 
tion (|36|) as in the XXZ chain. The important difference with XXZ is that there 
are now two kinds (0, 1) of quasi-particles, and their dressed scattering amplitudes 
depend on the type of quasi-particles. The BAE in the presence of holes read 



k' = 27r(p" + pI) -K°-kp° - A'-i * pi 
k' ^ 27r(pi + pI) - *p° - K"*pK 

We put them in "diagonal form" : 

2k' = 27r(p+ + p+) -K+*p+ 
= 27r{p- + pl)- K- -kp- . 

where p^ — p'^ ± p^. We can then express the BAE in terms of the p^ : 



2k' = 27^(p^ 
= 2TT{p- 



Ph_. 
Ph. 



where := {2tt - K^y^ -k {-2ttK^). Finally, we get the dressed BAE: 

k' ^2ir{p'^+pi)^ J2 ^"-'^pi, 



(76) 



(77) 



(78) 



(79) 



6=0,1 



where (f>° := i((f>+ + (f> ) and $i = $ ^ := i(<f>+ - <f> ). The scattering amplitude 
between particles of types a and b is then 

nA-A' 

5'^'''(A,A') = exp ' 



/•A — A 

Jo 



(80) 



We give the kernels $° in Fourier space: 

-27rcosh(7r — 37)cli ^^.^ 



$°(cj) = 



2sinh7w sinh(7r — 27)0; 



27rcosh7a; 



2 sinh 70; sinh(7r — 27)0; 



(81) 



5.5. CFT spectrum 

The CFT describing the staggered vertex model is a combination of two bosons (j>i 
and 02, with action (see dUl) 



-^-^/d'a^ [ffi(V<^i)'+52(V(/.2)'] 



2^, 



where 



9i 



go~^ 0. 



(82) 
(83) 
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As it is done for the XXZ model in Section 14731 the above action is inferred from 
the central charge and operator content of the BAE. Since the ground state is of 
the XXZ type and the Fermi velocity has the XXZ value, we easily find the central 
charge of the staggered vertex model: 



c = 2. 



(84) 



Similarly to the XXZ model, the primary operators correspond to higher magneti- 
zation sectors (magnetic operators) and global shifts of the Bethe integers (electric 
operators). In the present situation, each set of roots {A°},{Aj} can be excited 
independently, and we denote wPj-m} (resp. e°,e^) the magnetic (resp. electric) 
charges. The Wiener-Hopf technique^ can be adapte to the Z2 BAE, and we 
obtain the conformal dimensions 




^2 



1 




(85) 



(86) 



where 
kernels: 



± e^,m^ :— wl' ± m^, and gi^gi are encoded in the scattering 



1 

^^ = 4 



1 



j^+(0) 
27r 



1 



1 



27r 



The above spectrum corresponds exactly to two decoupled bosons as in 
ever, the electromagnetic charges obey the parity conditions 



+ e: 



[2], 



0[2] 



(87) 



. How- 



(88) 



which introduce a coupling of the excitation sectors (or boundary conditions) for the 
two bosons. The vanishing of 32 makes 02 a non-compact boson in the continuum 
limit, and each conformal exponent has a continuum of m" magnetic states above 
it. The level density of this continuum is given by the finite-size scaling of 32- The 
latter can be calculated numerically for large but finite N, by solving the BAE in 
sectors with 7^ 0. One gets the following results 



g2iN) cx 



for < 7 < f , 
for 7 ^ I , 



(89) 



where Nq is a 7-dependent cut-off value. The numerical results for 52 are shown in 
Fig.H 



5.6. Effective central charge and critical exponents 

In the presence of a twist ip, the conformal dimensions (|55|) - (I55]) are modified by 
the change e"*" — >■ e"^ -t- The ground state of the critical AF Potts model has a 



I 



I 
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twist tf = and we denote eo := ^- The effective central charge is thus 

Sen 6 , , 

Ceff = 2 ^ = 2 - - . (90) 

91 t 

Similarly to Section 14.41 some critical exponents of the Potts model at the AF 
transition are obtained from various sectors of the Z2 vertex model. 

• The £-\eg watermelon exponent has magnetic charge m+ = |. The anti- 
symmetric magnetic charge does not contribute to the conformal dimension, 
since 52 = 0. Hence: 

^ ^ g,{e/2y _ J_ ^ ^2_4 
^ 2 251 8t ■ ^ ^ 

• The 'thermal sector' with m+ = contains a state with twist (p = 7r(eo — 1), 
and exponent: 

^ 2gi 2 ^ ^ 

However, the first excited state in the Potts model has an exponent 
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which is not part of the electromagnetic spectrum (|85|) - ([85|) . This exponent 
may be interpreted from parafermionic CFTs for integer t (see Section [T]) , 
and is vahd for generic t. 
• In the case of the magnetic exponent Xh, the expected value when setting 
the twist to Lo — would be Xh — — ^ = I — ^ , but transfer-matrix 

1^ 2 " 2gi 2gi 8 2t ' 

calculations do not match this value. Rather, the magnetic exponent is 

Again, this value is consistent with parafermionic CFT for integer t. 

5.7. The limit Q 

We consider the following limit: 

TT 

7 = — + /i , u — hw , w fixed , /i — ^ . (95) 

In this limit, the double-edge _/?-matrix for the loop model has the form: 
TZ 

-7^ — > (1-w) 1 + w E + w(l-w) P := UBiw) , (96) 

where 1, E and P are defined by the diagrams in Fig. [5] The operators Ej and Pj 
generate a Braucr algebra with loop weight 0, defined by the relations 

P,P,±iP, ^ Pj±iP,P,±i (E,E,±,E,^E, 

P/ - 1 <^ E^ = (97) 

P,Pk = PkPj for |j - fc| > 1 [ EjEk = EuEj for |j - fc| > 1 

PjEj = EjPj = Ej 

EjPj±iPj = E,E,±, (98) 
Pj±iPjEj±i = E,E,±i 

Moreover, in the above limit, the vertex model is equivalent to the supersym- 
metric OSP(2|2) representation of the Brauer algebra. As a check, we can derive 
the Q ^ limit of the BAE, where A° hrj'^: 



-4 
"(< 



a=0,l j = l 

which are indeed the OSP(2|2) BAE. 

Loop models based on the Brauer algebra are integrable, but do not ad- 
mit a Co ulom b-gas description because loop intersections are present. In fact, it 
was showrPSI that they correspond to a universality class distinct from the 0{n) 
model. For integer loop weig ht, they admit a Bethe-Ansatz solution related to the 
0SP(n|2m) superalgebra 



I 
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P = 




Fig. 5. Generators of the Brauer model, built from the TL2jv('^ = 0) algebra. 



6. The anisotropic critical regime 

In this Section, we describe the critical regime of the Potts model based on the 
Z2 staggered 6V model (like the AF critical line), in the range < m < 7 where 
the coupling constants xi and X2 have opposite signJISl (see Tabled]). We refer to 
this regime as the anisotropic critical regime. Since it is based on the same vertex 
model as the AF critical line, the Bethe-Ansatz analysis is very similar, and so we 
will simply emphasize the differences with the previous Section. The resulting CFT 
is also a combination of two decoupled bosons (jji and 4>2, but in the present case, 
both bosons have a finite compactification radius. It is possible to write the toroidal 
partition sum as a generating series of conformal dimensions, and interpret 4)2 as a 
combination of two pairs of Majorana fermions. 



6.1. Bethe Ansatz Equations 

The Bethe roots are of two kinds: 

a° = A° , a] = \] + ZTT , 
and the BAE and energies read 



(101) 



exp [iNk(\';)\ = {-ly^-^e^^^ \{ [| exp [-ie°~^(A° - \\ 



6=0,1 1=1 



(102) 



where: 



^ 3> ^ cosh2A? - cos27 

a=0,lj = l a=Oaj = l 3 ' 



k{\) :— —i log 



sinh(i7 — A) 
sinh(i7 + A) ' 



e"(A) :=-*log 



sinh(»7 + I) o±iQV- ^ cosh(»7+|) 

— — ^, U {A}.- -I log — j- 

sinh(27 — ^) cosh(i7 — ^) 



(103) 

(104) 
(105) 
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States where {A"} = {A]} are solutions of an XXZ model with ^ = 27. The 
gromid state is a double Fermi sea with central charge c — 2. Quasi-particles are of 
two types 0, 1, with dispersion relation p6|) . and scattering kernels 



^ ^ -2vrsinh(7r-37)a. ^ 2™h7^^ 

2 cosh 70; sinh(7r — 27)0; ' 2 cosh 70; sinh(7r — 27)0; 



6.2. CFT spectrum and continuum partition function 

The conformal spectrum has the same form (|85p -(|86 |) as for the AF critical point: 

/^,^=i(^+m+V2^) > (107) 

= i - -+72^)% \ - rr^-^:j , (108) 

now with finite coupling constants 

The parity conditions also hold: 

e+ + = [2] , m+ + m" = [2] . (110) 

In this context, the toroidal partition function can be written as a generating func- 
tion for the conformal weights, and it reflects the coupling of the BC: 



1 



istag. 6V - ^;7;;Yk q"..^ 9''.-- . (Ill) 

e++e"=0 [2] 
m++m"=0 [2] 



The factors associated to charges e~ ,m~ are bosonic partition sums with coupling 
constant (72 = which can be written as products of two Ising partition sums 
with various BC. We introduce the Ising partition sum Z^^r' 1 with BC on the Ising 
spins (T — > (— 1)''(T, a — >■ (— 1)*^ cr, and the bosonic partition sum Zm,m'{g), with BC 
— > (/) + 27rm, — >■ + 211771' . The partition fvmction read^l^l 

^stag. 6V - ^ (-1)'^^''^+''^'-^ Zr,y^ Z^^y^ Z™,™'(.9l) • (112) 

m=ri+r2 [2] 
m'=r[+r'^ [2] 

In this form, the Z2 model appears as two Ising models and one compact boson, 
decoupled in the bulk and coupled through their BC. 
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6.3. Effective central charge and critical exponents 

In the presence of a twist (p, the conformal dimensions (|107|) - (|107|) are modified by 
the change e+ — > e+ + Like for the AF line, we denote cq :— ^, and the ground 
state has a twist (p = ttcq. The effective central charge is thus 

6e^ 12 
Ceff = 2 - ^ = 2 - — — . 113 
gi t(t - 2) 

• The £-leg watermelon exponent has magnetic charge m+ = |. Because of 
the parity condition (|110l) . the lowest possible value for the antisymmetric 
magnetic charge to^ is for even m+ and 1 for odd m"*". As a consequence, 
the watermelon exponents are 

\simf-± iff^o[4] , , 

• In the thermal sector with twist (p — cq, the energy exponent corresponds 
to charges e"*" = 2eo — 2, e~ — 0, which gives the constant value 

Xr=^^^^^^^l. (115) 

This is equivalent to the energy excitation in one of the Ising models. 

• The magnetic exponent Xh is given by the twist = ^. It turns out 
that the lowest exponent in this twist sector has electric charges = 
^ — 1, — 1. The exponent is thus 

(116) 



7. Conclusion 

The Q-state Potts model is a simple lattice model, in which the ferromagnetic 
critical point provides a basic realization of minimal CFTs. The antiferromagnetic 
part of its phase diagram has been less studied than the ferromagnetic one, but it 
contains three critical regimes of physical interest for Statistical Mechanics and Con- 
densed Matter Theory. In this review, we have shown that, through an equivalence 
to loop and vertex models, these critical points can be solved by Bethe Ansatz. This 
approach is powerful enough to discover the CFT spectrum, the S'-matrix descrip- 
tion of quasi-excitations, and also the critical exponents in the original statistical 
model. 

Most of the material exposed in this review is valid for generic values of Q. When 
Q is a Beraha number {Q — A cos^ j , with t integer) , the Potts model can be mapped 
to an RSOS height model (see Section Although the partition functions of the 
two models are identical (even at the lattice level), the excitation spectra can be 
very different. The essential mechanism at play is that some excited states of the 
vertex model have a vanishing contribution to the RSOS partition function. In 
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the CFT language, the ehmination of these states amounts to null-state equations, 
which are the basis for the determination of minimal series of CFTs. RSOS models 
with homogeneous Boltzmann weights (corresponding to the self-dual Potts model) 
were shown to have four different physical phases I, II, III, IV, with critical transitions 
separating phases III-IV and I-II. These two transitions are lattice realizations of 
the Minimal and Z^-parafermionic (for k ~ t — 2) CFTs, respectivelj^^. 

Back to the Potts model with Beraha values, this tells us that the RSOS version 
of the ferromagnetic critical point relates to a Minimal CFT, whereas that of the 
BK phase relates to the Z^.-parafermionic CFT. The question of the RSOS version 
for the AF critical line is more difficult, but it has been found recentljKI that this 
is also described by the Z^-parafermionic CFT. This relation is supported by the 
study of critical exponents and of the algebra of chiral currents. 

Although the critical phases of the antiferromagnetic Potts model are now iden- 
tified and well studied, there remain interesting open questions to the mathematical 
physicist. 

From the point of view of intcgrable models, the understanding of the Z2 model 
would be more complete if one could extend the relation to braiding algebras (see 
Section [577| . Also, an analytic derivation of the effective coupling constant for the 
non-compact boson should be possible. Moreover, one can generalize the Z2 con- 
struction to a Zp six-vertex model, where the Zp symmetry will be reflected in the 
whole Bethe-Ansatz structure. This opens a new series of models, similarly to what 
appears in the fusion procedure. 

Finally, the Potts model is a good framework to study the physics of cluster 
or spin interfaces in lattice models. Current studies focus on the links between 
boundary CFT, Schramm-Lowner Evolution (SLE) and free massless field theories. 
In the case of the antiferromagnetic Potts model, the Fortuin-Kasteleyn cluster 
boundaries are well described by the ^-leg watermelon exponents, all related to the 
first boson 0i, but the geometric expression of the second boson (/)2 is lacking. This 
additional degree of freed om ma y lead to an extended version of the SLE process, 
as was introduced recentlji^^^^. 
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